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ABSTRACT 


In order to simulate the space environment for basic research into the crystal growth mechanism, Hg 0 .sCd 0 . 2 Te 
crystals were grown by the vertical Bridgman-Stockbarger method in the presence of an applied axial magnetic 
field. The influence of convection, by magneto hydrodynamic damping, on mass transfer in the melt and segre- 
gation at the solid-liquid interface was investigated by measuring the axial and radial compositional variations 
in the grown samples. The reduction of convective mixing in the melt through the application of the magnetic 
field is found to have a large effect on radial segregation and interface morphology in the grown crystals. Direct 
comparisons are made with a Hg 0 .sCdc. 2 Te crystal grown without field and also in the microgravity environment 
of space during the second United States Microgravity Payload Mission (USMP-2). 


1 Introduction 


The Hg based Il-VI semiconductor, Hgi- x Cd x Te , is important for application of infrared detection and 
imaging applications. Its band gap, Eg , varies linearly with composition, x, 

E g = 1.9x - 0.3 (!) 


which allows it to be compositionally tuned for a broad range of wavelengths from 0.8 pm to the far-infrared 
spectrum beyond 30 pm. Due to the high Hg vapor pressure present during the processing of this material it 
must be sealed in a thick silica ampoule. This causes a large destabilizing horizontal temperature gradient to be 
generated due to a difference in the thermal conductivities of the melt and crystal and the release of latent heat 
at the growth interface in the presence of the silica ampoules as predicted by Jasinski and Witt. 1 Additionally, a 
stabilizing vertical solutal gradient is produced by the rejection of the denser solute, HgTe, in to the melt. These 
phenomena, coupled with a large solutal-to-thermal expansion coefficient 
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( 3 ) 


and a large thermal-to-solutal diffusion coefficient ratio 


give rise to double diffusive convection during melt growth. This buoyant convection, due to fluid density dif- 
ferences which are often small, can produce strong convection in all noncapillary solidification processes. The 
equation which determines the intensity of convection is the Navier-Stokes (N-S) or momentum equation, 


p^ = -VP + pV 2 V + pG. (4) 

Expanding equation 4 in the x direction with the Boussinesq approximation for the density dependence on 
temperature and solute concentration yields, 
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where p is the fluid density and p is the fluid viscosity, g is the acceleration due to gravity, 3 t is the thermal 
expansion coefficient, 8 S is the solutal expansion coefficient, AT is the driving temperature difference and AC is 
the driving solutal difference. This equation is nothing more than Newtons first law: 


Dt ^ 


(6) 


i.e., the summation of forces are equal to the product of an objects mass, Af, and acceleration, ^ , derived for 
an elemental fluid volume. Therefore, every term in the N-S equation is merely a force. The driving temperature 
difference, AT, in the thermal buoyant force is determined by the energy equation: 


= KV2T 


(7) 


where p is the density, C p is the specific heat at constant pressure and K is the thermal conductivity. Similarly, 
the driving solutal difference, AC, in the solutal force is determined by the species equation: 
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( 8 ) 


where D is the solutal diffusivity. 


During crystal growth, the severity of convection is governed by the highly coupled form of equations 4, 7 and 
8. Namely, ratio of the thermal and solutal buoyant forces to the viscous forces. The ratio of these forces is given 
by a dimensionless parameter known as the Grashof number, Gr, 


Gr t = 


gftrATI 3 


and Gr , 
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(9) 


where L 3 is the characteristic length scale of the system (for vertical Bridgman-Stockbarger growth the length 
scale is the crystal radius, r cr ), v is the kinematic viscosity and a is the thermal diffusivity. The subscripts, t and 
s, denote thermal and solutal respectively. 


As mentioned previously, the thermal force is due to density differences arising from unavoidable radial tem- 
perature gradients in the melt which promote convection. Solutal forces are due to density differences due to the 
rejection of solute during solidification. In dilute alloys the solute rejection has little influence on density and 
the Gr s is small compared to the Gr t . However, in nondilute alloys such as Hg\- x Cd z Te this solutal force is 
significant and has a large effect on the fluid flow intensity and fluid flow pattern in the melt which can, in turn, 
have a significant effect on crystal composition. This is easily seen from equation 1. A few percent change in 
composition, x, will dramatically change the samples bandgap, E g . 
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1.1 Methods to Reduce Segregation 


It has been shown that fluid flow can interfere with segregation and lead to the nonuniform incorporation 
of conservative and nonconservative dopants into the grown crystal. If the dopant distribution is not radially 
uniform within a few percent the material is useless to device manufacturers. In order to control convection 
and segregation in the material, the growth system parameters such as charge size, temperature gradients and 
acceleration due to gravity must be modified. Through, for example, growing small diameter crystals, processing 
in space or under the influence of applied magnetic fields. 


1.1.1 Reduction of Charge Diameter 


As the severity of convection in crystal growth from the melt is controlled by the Grashof number shown in 
equation 9 which has a cubic dependence on the crystal radius, r rc , the radius is one of the obvious parameters 
to change in order to reduce convection. To this end, experiments were performed with various charge diameters 
during vertical Bridgman growth of Ga-doped Ge. 2 Results showed that the crystal diameter had to be reduced to 
capillary size to avoid convection in the melt. Obviously, capillary sized crystals are of no value to the electronic 
industry. They have, however, been found useful for important studies of the semiconductors thermophysical 
properties such as diffusion coefficient measurements. 


1.1.2 Reduction of Gravity 


Another parameter which can be used to control convection is the acceleration due to gravity, g . Experi- 
ments 3 — 8 performed in the reduced gravity environment of space have shown that diffusion limited growth can 
be obtained. The real advantage of processing in microgravity is that it provides the opportunity to grow larger 
diameter crystals in the diffusion limited regime than is possible on earth. This can be readily shown through 
scaling analysis by comparing the convective acceleration force to the thermal driving force. This relationship is 
shown in figure la as the Reynolds number, Re, 


Re — 



( 10 ) 


verses the Grashof number, 
the Schmidt number, Sc = 
( ReSc ) = Pe s . Noting that 


Gr t . Multiplying the Re by the Prandtl number, Pr — — , for the first case, and by 
for the second case, yields plots similar to figure la where ( RePr ) = Pet and 
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( 11 ) 


for typical semiconductors, a Re of 10 in figure la becomes a Pe t of 0.1 and a Pe $ of 100 in figures lb and lc 
respectively. 


These figures show that diffusion limited heat transfer ( Pet < 1) can be obtained at high values of g (say on 
earth) and at the same g - level have convection dominated mass transfer (Pe s > 1). Therefore, to get diffusion 
limited mass transfer the Gr t (g-level) has to be reduced to very low levels such that Pe s < 1, at which point 
Pet is very small (recall* the ratio of the two was 10 3 ). It can be shown that the functional dependence of Pe s 
on Gr t is 

Vr v 

Pe s = — — « fen 
v D 

The quantities 5,AT,z/, a are fixed for a given crystal growth system and on earth g is fixed which leaves the 
crystal radius as the only adjustable parameter. However, in space one can lower g easily by four orders of 
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Figure 1: (a) Plot of Re vs. Gr t , (b) Plot of Pet vs. Gr i? (c) Plot of Pe s vs. Gr t for a typical semiconductor 
material. 


magnitude. Thus, from equation 12 (roughly speaking), with all other parameters constant 

r space ^ t 9space 
T earth \9earth 

i.e., space provides an excellent opportunity to increase the yield of the grown crystals. 


(10)3 « 30 


(13) 


1.1.3 Application of Static Magnetic Fields 


A less expensive v/ay to modify convection is through the application of magnetic fields. The basic mechanism 
for the interaction of an applied magnetic field and a molten semiconductor involves the electrical currents induced 
by the movement of a conductor in the presence of a magnetic field. When convection is present, flow of this 
electrically conductive fluid perpendicular to the field results in a Lorentz force that opposes that motion. To 
determine the effectiveness of the magnetic field on the suppression of convection it is necessary to solve Maxwell’s 
equations and the fluid flow equations for the velocity and electromagnetic fields. A complete derivation of the 
MHD equations is given in 9 and basically results in the incorporation of the Lorentz force 


F L orentz=v(E + V XB)XB (14) 

into the Navier-Stokes equation (eq. 4). The neglection of the induced EMF E has been shown to be valid in the 
limit of two-dimensional flow. 10 With this approximation equation 14 reduces to 


F Lorentz = <t{V X B) X B (15) 

where a is the electrical conductivity of the fluid, V is the fluid velocity vector and B determines the magnitude 
and orientation of the magnetic field. In many respects the damping of the convection by magnetic fields is similar 
to a change in the viscosity of the liquid- This has given rise to the term ’magnetic viscosity 5 which is analogous 
to the kinematic viscosity and is commonly denoted as 

Vmagn ~ 0 (16) 


where //o is the fluid permeability. 
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The nondimensional parameter which determines the effectiveness of the magnetic field on the suppression of 
convection is the Hartmann number, Ha, 


Ha = BL 


(17) 


which is the ratio of the Lorentz force to the viscous force. For typical semiconductors and medium to large field 
strengths, Ha as 10 6 , indicating that the electrical currents are better at transporting energy over larger distances 
in the melt than the viscous forces. Using a similar scaling argument as that presented above and noting, 


Vmagn 

V 


10 6 


(18) 


it can be shown that 



(19) 


Thus, both processing in microgravity and with the application of a magnetic field provide the potential to 
increase the quality and yield of semiconductor material. 


This research focuses on the latter approach to suppress convection in a molten semiconductor. In the next 
section, a brief description of the experimental magnetic crystal growth system will be presented. Then we will 
discuss crystal growth results in a magnetic field. Direct comparisons will be made between microgravity and 
zero magnetic field cases. 


2 Crystal Growth Apparatus 


The ingots were grown by directional solidification in a vertical Bridgman-Stockbarger furnace. The furnace 
consisted of: a hot zone where the sample was melted; a booster heater located between the hot zone and gradient 
zone and a 0.3 cm heat extraction plate placed between the gradient and cold zone, which produced the high 
axial gradients necessary for this research. A guard heater was utilized on the end of the hot and cold zones 
to minimize any thermal losses out the end of the furnace. The furnace surface was maintained at 20 C with 
a water coolant loop. This furnace was designed to provide similar thermal environments as those produced by- 
NASA’s two microgravity furnaces, the Advanced Automated Directional Solidification Furnace (AADSF) and 
the Crystal Growth Furnace (CGF). The ampoule was supported by a cartridge assembly described in detail 
in. 11 This configuration was used during microgravity growth of Hgo.sCdo.-iTe and was adapted for this research 
in order to have a direct comparison. Crystal growth was achieved by translating the furnace which allowed 
crystall iz a tion to take place in the homogenous region of the magnetic field. The furnace assembly was coaxially 
aligned with a superconducting magnet system designed specifically for this research 12 as shown in figure 2. 


3 Experimental results and discussion 


The effect of convection on segregation was determined by measuring radial and axial compositional variations. 
Before sectioning, surface morphology and axial compositional measurements were made on the as grown ingots. 
This data, along with growth orientation, was used to determine the sectioning of the ingot for further study. 

For radial segregation measurements wafers were cut perpendicular to the growth direction. Axial segregation 
was determined from centerline samples cut parallel with the growth direction. The macrosegregation of CdTe 
was determined by Precision Density Measurements (PDM) 13 while microsegregation was determined from the 
transmission edge of the IR transmission spectrum and Energy Dispersive X-ray spectroscopy analysis (EDX). 
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Figure 2: Schematic of superconducting magnet and crystal growth system. 


3.1 Axial compositional distribution 

For all growths with zero magnetic field the axial compositional profiles compared well with the profile pre- 
dicted by the one-dimensional model by Clayton and colleagues 10 see figure 3. The axial compositional profiles for 
all growths under the influence of a magnetic field also compared well with the one-dimensional diffusion model 
with the exception of an abrupt compositional rise near the end of the ingot. This rise in composition has been 
attributed to constitutional supercooling and has been described in detail in. 11 



0.0 2.0 4.0 6.0 8.0 10.0 12.0 14.0 


Axial Distance (cm) 

Figure 3: Comparison of experimental axial compositional distributions to one dimensional theory. 
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3.2 Radial compositional distribution 


Radial compositional profiles for growth in 0, 5, 10 and 50 k G magnetic fields are shown in figure 4. All 
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Figure 4: Radial compositional distributions. 

wafers were taken at an axial position of 60 mm from the first to freeze portion of the as grown ingot. The 
pattern of low CdTe content at the center of the wafer for zero magnetic field is consistent with other results 
obtained previously in the absence of a magnetic field. This large radial compositional variation is consistent with 
a concave (toward the solid) interface curvature w’here the less dense component (HgTe) tends to accumulate at 
the lowest gravitational potential. There is a remarkable improvement in radial segregation in all samples grown 
in the magnetic field which indicates that reduction of convection through the application of a magnetic field is 
advantageous for maintaining the solid-liquid interface shape required to minimize the compositional variation 
transverse to the crystal growth direction. 

To further investigate the effects of the magnetic field on interface shape, experiments were performed to mark 
the solid-liquid interface during growth with and without the influence of the stabilizing magnetic field. This was 
accomplished by quenching the sample after approximately 65 mm of growth. The resulting interfaces are shown 
in figure 5. The interface curvature for the crystal grown in a 50 kG field was a factor of three less them the 
interface obtained with zero field. This suggests that the convection has a significant influence on the interface 
morphology HgCdTe system. Information on the fluid flow velocities and fluid flow patterns are necessary to 
make definitive conclusions about the convective effects on interface morphology. 

To estimate the relative magnitude of the convection in the HgCdTe system the analytical model for double- 
diffusive convection by Hart 16 was utilized. This model was originally developed to study salinity gradients 
but has been shown to be applicable to nondilute crystal growth systems. 17 Hart solves the linearized forms 
of equations 4, 7 and 8 in terms of a single solution parameter dependent only on the solutal Raleigh number, 
M = (-0.25 Ra $ )*. Using this analysis and the parameters for the HgCdTe system, the maximum velocity in the 
system without the magnetic field is 10 /zm/sec. This is 5 times larger than the velocity required for diffusion 
limited growth (2 jzm/sec) as determined from a Pe s 1 calculation and results in large convective interference 
with segregation. The maximum velocity during crystal growth with magnetic field is on the same order as the 
diffusion limited velocity and corresponding radial segregation is nearly uniform. It is of interest to note that if 
the flows in the HgCdTe system are indeed as slow as predicted by this analysis a smaller magnetic field should 
still produce uniform radial segregation. To this end, growths with magnetic field strengths of 10 and 5 kG 
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Figure 5: Solid-liquid interface shapes for growth (a) zero magnetic field, (b) 50 kG magnetic field. 


were performed. This represents an order of magnitude change in field strength and in all cases uniform radial 
segregation was obtained as shown in figure 4. 

To further quantify the effects of the magnetic field on the radial segregation, the experimental results were 
compared to the analysis of the lateral solute segregation associated with a curved solid-liquid interface during 
steady-state unidirectional solidification of a binary alloy as derived by Coriell and Sekerka. Briefly, they solved 
the species continuity equation for the radial solute concentration in the solid crystal at the solid-liquid interface 
with the assumptions of no convection in the liquid and that the solid-liquid interface could be represented by 
a Fourier series. Their work showed that in the limiting case (i.e., when the solute boundary layer is thick 
compared to the deviation of the solid- liquid interface from planarity), the transverse segregation in the solid 
is proportional to the deviation of the interface from planarity, the proportionally factor being just the product 
of the unperturbed concentration gradient and the distribution coefficient. The Coriell formalism is developed 
for a 2-D rectangular geometry which does not completely describe the crystal growth geometry in the HgCdTe 
system. To rectify this problem a numerical code based on the Coriell formalism was developed to predict radial 
compositional distributions in HgCdTe. The configuration used was a cylindrical system at steady state growth 
conditions and the solid-liquid interface was expressed in terms of a series combination of Bessel’s functions. The 
details of this analysis are described in detail elsewhere. 19 These calculations are applicable for crystal growth in 
a stabilizing magnetic field where the field strength is sufficient to suppress convection. The radial segregation 
for growth in a 50 kG magnetic field is shown in figure 6. There is excellent agreement with the Coriell formalism 
and indicates that the radial segregation is approaching the diffusion limited regime. 


4 Comparison with microgravity crystal growth results 


A Hg 0 .gCd 0 .- 2 Te alloy crystal was grown in the microgravity environment aboard the space shuttle during 
the Second United States Microgravity Mission in March 1994. The grown ingot was approximately 16 cm long 
and 0.72 cm in diameter. During the mission portions of ingot were grown while the orbiter was in different 
attitudes. The influence of orbiter attitude and, hence, the residual drag acceleration on the growth process was 
investigated by examining the radial composition of the ingot during the +YVV/-ZLV and -ZVV/-XLV portions 
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Figure 6: Comparison of radial compositional profile obtained with the application of a 50 kG magnetic field to 
theory. 


of the mission see figure 7. 

The radial compositional profiles for the different attitudes are shown in figure 8. These measurements were 
taken at 15 degree intervals at circumferences within the sample at progressively increasing distances from the 
edge, and starting at 100 fj.m in. As can be seen, the asymmetry in the compositional variation across the 63 mm 
wafer shows a definite correlation with the direction of the residual acceleration vector resolved on this plane. This 
is attributed to residual fluid flow. This growth configuration is unstable in that there is a residual acceleration 
component in the direction of solid to liquid, which is destabilizing. In similar vein, the sample cut at 106 mm also 
shows asymmetry, and once again the effect of the residual acceleration component is to affect the distribution 
of the liquid such that the heavier material is sedimented asymmetrically by the residual force. In this part of 
the growth, there was a stabilizing vector, albeit small, in the direction of liquid to solid. This contrasts with 
the 63 mm wafer, and has resulted in a more homogeneous wafer. The compositional profile for growth in a 50 
kG magnetic field is shown in figure 8c. The compositional distribution is more uniform and symmetric which 
indicates that the magnetic field is efficacious in suppressing convection that may result from any misalignment 
with the gravity vector under Ig conditions. 


5 Summary 


To summarize, the most interesting finding of this research is that convective interference with segregation 
and the interface morphology in the HgCdTe system is quite subtle. By reducing fluid flow from 10)um/sec 
(no magnetic field) to 2^m/sec (magnetic field), the solid-liquid interface deviation from planarity was reduced 
by a factor of 3. This is attributed to the large separation of the liquidus and solidus phases in the HgCdTe 
system which produces a highly compositional dependent solid-liquid interface. The importance of the solutal 
term compared to the thermal term on density should be noted. A change of 5 to 10 percent in composition 
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Figure 7: Components of residual acceleration at AADSF during USMP-2. 


has the same effect on density as a 150 C change in temperature. Therefore, a small reduction in fluid flow 
(i.e., reducing the build-up of the solute at the sample center) will tend to produce a more planer solid-liquid 
interface. Obviously, one can never achieve a planer interface due to the differences of thermal conductivities of 
the melt, solid and containing crucible. However, this effect can be minimized by judiciously choosing the thermal 
boundary conditions. 


In microgravity crystal growth a strong correlation between the relative magnitude and orientation of the 
residual acceleration due to drag was observed. If the drag vector is not aligned parallel to the growth direction 
(liquid to solid) three dimensional fluid flow will result which destroys the symmetry of the system and leads to 
nonuniform radial segregation. 
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Figure 8: (a) Radial compositional map for the +YVV/-ZLV attitude (63 mm), (b) radial compositional map for 
the -ZVV/-XLV attitude (106 mm), (c) radial compositional map for a sample grown in a 50 kG magnetic field 
(69 mm). 
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